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ABSTRACT Many molecutar motort movo unkllracttonaily along a DNA strand powered by nucleotide hydrocele, Thooo 
rootore art muKJroortc ATP at at with mora than ona hydrocele efte. Wo prat ant hara a modal for how thaaa motort generate 
tha raqulatta foroa to prooaaa along thalr DNA trade TWt novel maohanlam for foroa generation la baaad on a fluctuating 
al actroatattc INMd d ^^^ udi0 ^^'^^ Wa apply tha prlnolpla to explain tha motion of oartain DNA hailoaaat and 
tha portal protein, tha motor that bactarlophagaa uaa to pump tha ganoma Into thalr oa p a lda. Tha motor can ravaraa Its 
a (racoon wnnouc ravaraing ma potarny ot rtt atacTroataiic new, mat la, wttnout major structural moomoationa of ma protam. 
Wa alao ahow that tha motor oan b« drlvan by an lm 
motor and to ATP aynthaaa. 



INTRODUCTION 

A variety of protein motort have the property that they 
process unidirectionally along a polymer track, driven by 
nucleotide hydrolysis. Actin and tubulin provide the 
tracks for the moat well-studied examples, myosin, kine- 
sin, and dynein. DNA is also traversed by a medley of 
protein motors that replicate, repair, and package the 
genome. These protein machines differ from myosin and 
kinesin in several ways that suggest their mechanical 
mechanism is different. First, certain DNA motors appear 
to be rotary, rather than reciprocating; i.e., they process 
not by "walking" along their track, but by rotating about 
it, so that their procession is helical rather than linear. 
Moreover, they are multimeric and may hydrolyze nucle- 
otides at several catalytic sites (Patel and Hingorani, 
1994; Patel et al., 1994). Although it is clear that these 
motors are driven by nucleotide hydrolysis, the precise 
mechanism by which phosphate bond energy is trans- 
duced into a directed force remains mysterious. Here we 
propose a novel mechanism for rotary molecular motors 
that we apply to a model for the portal protein and certain 
DNA helicases. 

For many bacteriophages, the culmination of assembly 
involves the insertion of several hundred kiiobases of 
genomic DNA into the capsid head. This is an astounding 
packaging job, for the length of the DNA is -3500 times 
as long as the capsid, yet the entire process is completed 
within 15-90 s, which suggests a velocity of more than 1 
pm/s. Although the packaging mechanism remains 
largely mysterious; several known features provide some 
clues. First, one ATP is hydrolyzed in the packaging of 
two base pairs, or 9 X 10 3 ATPs are needed to package 



the 1.8 x 10 4 base pairs in the phage $29 (Ouo et al., 
1987). This suggests that DNA is driven into the capsid 
by some sort of ATPase motor. Second, DNA is inserted 
through a protein complex at the capsid base, the portal 
protein, or -connector" (Black, 1988; Barnshaw and 
Casjens, 1980; Valpuesta and Carrascosa, 1994). The 
molecular geometry of this structure is known for several 
phage species (Carazo et al., 1986; Casjens et.al., 1992; 
Donate and Carrascosa, 1991; Dube et al., 1993; 
Valpuesta et al., 1992), and its geometry and symmetries 
suggest that DNA may be inserted in a rotary and/or 
transitional motion through the central channel. 

DNA helicases are also nucleotide-driven motors that 
can translocate along double- or single-stranded DNA at 
speeds of up to 1000 nucleotides/s ~ 34 nm/s. Although 
rolling or walking mechanisms have been suggested for 
certain helicases (Lohman, 1992; Wong and Lohman, 
1992), die structure of some helicases is similar to, but 
simpler than, that of the portal proteins, suggesting a 
rotary motion. For example, the T7 gene4 helicase is a 
doughnut-shaped hexameric ATPase with three hydroly- 
sis sites that translocates along ss-DNA (Egelman et al., 
1995; Hingorani and Patel, 1993; Patel and Hingorani, 
1994; Patel et al., 1994). Fig. 1 shows the geometry of a 
portal protein from T3 bacteriophage and the helicase 
from T7 bacteriophage. 

Here we shall present a model for the mechanochemistry 
of these two DMA motors, which we propose operate by 
similar mechanisms. For simplicity, we shall focus on the 
T7 helicase. 
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THE FLASHING FIELD MODEL 

The model is based on the following central assumption (see 
Fig. 1): Binding of ATP to a hydrolysis site induces a 
conformational change that exposes a pair of negative 
and positive charge regions near the inner surface of the 
channel 
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The chatged regions are not of equal size, and they are 
oriented at an angle to die ciicumferential meridian. We will 
assume that there is one charge pair region per nucleotide 
hydrolysis site, and that the axes of the charge pair regions are 
tilted with respect to the axis of the hole. The helicase binds 
three nucleotides at a time, and so we inter it has three hydro- 
lysis sites (Ptfel et aL, 1994; Patel and Hingorani, 1995); 
however, the portal protein has perhaps 12 or 13 hydrolysis 
sites, located m two tiers (Casjens et aL, 1992; Dube et aL, 
1993). As we will see, the pattern of hydrolysis affects the 
performance of the motor (analogous to the firing sequence of 
a car engine). Therefore, if each hydrolysis site controls the 
expression of a single charge pair then, tor more than a tew 
hydrolysis sites, the number of firing sequences becomes too 
large to compute. 

Qualitatively, the model works like this. Hie negatively 
charged phosphates spaced along the backbone of the DNA 
sequentially with the field of the charge pair 



regions as they "flash" on and off with the hmdi^g and hydro- 
lysis of ATP. Each charge pair field gives the closest phos- 
phates) an "electrostatic push** in the direction of the charge 
pair axis. The net effect of the charge pairs flashing on and off 
creates a sustained torsional and axial thrust It turns out that 
the order in which the fields flash is not too imp o rtant; even 
random flashing win drive the rotation. Note that the motor is 
an doUmpiutic machine that transduces strain flurtiiattons gen- 
erated by nucleotide hydrolysis into a fluctuating ekctric field 
Biased diffusion plays no role in torque generation; Brownian 
motion serves only as a *1ubricanT to ensure the rotor does not 
hang up in local energy minima. Fig. 2 illustrates the model 
geometry and the operating principle. 

APPROXIMATIONS AND PARAMETER VALUES 

In deriving the model equations we make several simplify- 
ing approximations. None of these are essential to the mech- 



2266 



Btorrvtaoai Joumri vttmtt Dmmbar 1006 




F10URE 2 (a) The potential (w&i /in*) and the force (dajAed /i**) of t linage off-«xii, screened dipote (Eq. 2) teen by a stof> negirive phosphate charje 
along the DNA becktaste. The dashed line U the force on a negative chart©, Id tbe shaded regions the force on the phosphate is to the right Althoogh there 
to • small regie* (whtit) where the force acts strongly to the left, the flashing, asymmetric potentisl produces a net farce to the right (n) Ittosttating the 
flashing field prindpto; A pair of nes^ve charges sre rig^y connected snd 

noise and dipote 1^2, or 3. The dipotes flash between coofiguratioos 1, X and 3 diner ■eqoeotiaBy or randomly, one being "on" at all times- The potentials 
are each repeated! periodically in the x direction. The flashing fields drive the negative charges to the right (c) The heticase baa been slit and onroOed to 
show the inner surface of the hole. The diagonal line is the phosphate b ackb o n e of the DNA; the p ho s phates are shown as Mack dots (•). In mis figure, 
the bote contains one tarn of the helix. The three dipotes, labeled A, B, and C, are spaced along the top of the bole and inclined with respect to the axis 
of the bote. Assume that dipote A flashes on; the backbone charge closest to A feels the dipolc field and is given both an axial and a cuuuufa e u nal posh. 
This brings another dprge into range of dipote B. If dipote B switches on after A switches off, this phosphate is given another push. Thus as the dipotes 
flash on they exerf both an angular and a longitudinal push on the doaest phosphates, bringing another phosphate within range of the next dipote's field. 
Because rc<atk»dfdiffiMk)o is very rapid, "dead spots" do not hang up the DNA motion; i*^ Browdan motion acts as a 'lubricant'* The uddirectionaliry 
is determined by'tJy o rien t ati on of the dipote Adds. The ability of the motor to reverse its direction of rotation requires an asymmetric field, which we m od el 
as a mritipote fieltfc^ 

positive charge is twice* the magnitude of the negative charge. The arrows depict the force experienced by a negative charge in the field. The potential profile 
in a corresponds so the section shown by the arrows. 



anism, but they arc necessary to formulating a tractable 
mathematical model. 

1. The DNA is stiff. This is a good approximation be- 
cause the perlisUnce length of DNA is about 50 nm, 
whereas the length of the DNA motors is £12 nm. 

2. The DNtil strand fits snugly into the hole in the motor. 
The diameter of a B-DNA strand is —1.8 nm, whereas the 
mameter of the Bole is 2.5 nm (belicase)-3.7 nm (portal 
protein); thus there may be considerable wobble accompa- 
nying the motion. This greatly increases the complexity of 



the computations without adding much to the principle of 
operation. Therefore, we have simplified the model by 
neglecting the wobble, e.g., the DNA threads through lite 
hole snugly. In fact, the effective diameter of B-DNA, 
including hydration and charge, may be larger than that 
measured diameter, and so the DNA strand may indeed fit 
snugly into the central passage. 

3. The sizes of the positive and negative charge regions 
in the charge pairs are unknown, and so for simplicity we 
have assumed a simple asymmetrical field consisting of a 
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dipole plot a monopole. (The recently elucidated struc- 
tint of the PI portion of ATP synthase shows a triplet of 
equally spaced, unequal, charged regions lining the ori- 
fice facing the 7 theft. Cf. figure 3 of Abrahams at al., 
1994.) Beeanae the angle of die field la unknown, we 
have assumed it la 45*. 

4. The die l e ctri c environment of the bole is unknown. 
Therefore, we have assumed a dielectric constant typical of 
protein interiors (• - 4-7). Hie central cavity in Fl-ATP 
synthase is hydrophobic (Abrahams et at, 1994), and we 
have taken this as justification for this assumption. 

5. Hie number of phosphates interacting with each 
charge pair depends on the fit of the DNA in the motor 
channel and the interior dielectric environment We have 
assumed that only two phosphates sense the electrostatic 
field at a time. Increasing the number of phosphate charges 
that sense the field will not alter the qualitative behavior of 
the model. 

6. The net hydrolysis rate is chosen arbitrarily at 
100/s. The total flashing rate for 429 portal protein can 
be obtained from the measured nucleotide hydrolysis 
rate, -1500 a" 1 . However, the number of ATPase sites is 
unknown. For helicase, the net hydrolysis rate has not yet 
been determined. However, knowing that the velocity is 
about 1000 nucleotides/s (-0.34 jim/s) and assuming 
completely tight coupling (i.e., one hydrolysis per nucle- 
otide), an upper bound on helicase hydrolysis can be 



are known, and the above 
in only demonstrate the prin- 
the model points to definite 



Until the detailed 
quantities are available, 
dple of operati o n; 
measurable properties. 



MODEL EQUATIONS 



the pair of Langevin equations: 



Torque due Load Brownian 
to potentials 



The equations describing the motion of the motor can be 
written in two equivalent ways: as stochastic ordinary dif- 
ferential (i.e., Langevin) equations describing the motion of 
a single motor, or as diffusion (i.e., Fokker-Planck) equa- 
tions describing die motion of an ensemble of motors 
(Doering, 1990; Gardiner, 1985; Riskin, 1989). The 
former are generally easier to understand, simulate nu- 
merically, and compare with experiments, whereas the 
latter are easier to obtain statistical properties of trajec- 
tories. We shall present the model as Langevin equations; 
the equivalent Fokker-Planck equations are presented in 
the Appendix, along with the computational algorithms we 
have employed. 

The equations describing the motor are as follows. Let 6 
be the cyclical coordinate of the position of a dipole on the 
inner surface of the helicase; We put our coordinate system 
on the DNA, so Chat it remains stationary while the helicase 
spins and translates. The circumferential ($) and transla- 
tions! (z) co m ponen t s of the helicase motion are given by 
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0*6*2* 



^dt — * — T ' + W* T -Uf)> 



Axial thrust Load Brownian 
due to 



(lb) 



Here die subscripts $ and 1 refer to rotational and longitu- 
dinal quantities, respectively. U/f) are gauasian white 
noise processes satisfying <«0) - 0, <fif) «f - i)> 
« 6Xf - s). The motor is driven by the time variation of the 
potentials' amplitudes, which are driven by the hydrolysis 
cycle. Although the magnitudes of the charge regions are 
not equal, for simplicity we shall begin by modeling the 
spatial dependence by a periodic array of dholes. Later we 
shall show the effect of the charge 1 



The dynamics of the rotational motion as described by Eqs. 
la and lb reveals the stochastic character of the motion, 
including forward and backward fluctuations. Fig. 3 a 
shows a stochastic simulation of the trajectory, (0(0, 40)* 
for the helicase (three-dipok) situation showing the helical 
progression. The stochastic simulation method is described 
in Appendix A and the diffusion equation method is de- 
scribed in Appendix B. 

The mechanical performance of the rotary motor is en- 
capsulated in its load-velocity curve (Svoboda and Block, 
1994), which is best computed from the Fokker-Planck 
representation of the dynamics. However, this c omputatio n 
is difficult if the axial and rotational motions are indepen- 
dent, and so we shall hereafter assume mat the helicase 
processes such mat the rotational and translations! motions 
are coupled; that is, the protein winds along the DNA like a 
not on a screw. This reduces the degrees of freedom to ooe, 
so that we can compute die rotary motion indepen de ntly , 
men obtain die translations! motion by tri gonom etry. The 
stochastic simulation in Fig. 3 a shows that the motor works 
without this simplification. 

For tins calculation we use die potential due to an "offnuas 
screened dipole potential" (Le~, the difference between two 
equal but oppositely charged exponentially screened Coulomb 
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the radius along the inner surface. The potential due to each 
dipolc interacting with the array of M charge* 00 the DNA 



V^)-iv^-2ir5r). (2b) 

The time-dependent potential VU 0 acting in Eq. 1 la, at 
each instant of time* that due to one or another of the N 
phase ■shifted dipoles 00 the inner surface of the hdicaaei 
The configuration of dipoles switches in time between the N 
phase-thifted configurations so that the full pote nt ial in the 
nth configuration can be written in terms of V WTty (0) ac- 
cording to 
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FIGURE 3 («) Sapchaatk trajectory computed from Bqs. 1-3 for the 
helkate with thrm flashing dipoles Intcnctinf with two ph o spris aos el a 
time. Dark Una - */) [redian]; gray line - 40 [mn). <*) Loed4onfW 
angular velocity cwrva computed from the dtffutioo oquationi described in 
Appendix B. Some simnlstkwis were alto canied out using a shooting 
method with the program xpp (available via WWW at http*y/info.pttLeda/ 
-phase/). Panwatera; radios - 8 nm; wefl depth - 8 *«r, twitching 
me - 100ft; tame screening length * SO; off-axis - OS; D $ - 7.5 x 
10* nm™/*.] (C at eadmated from the expression for the frictknal drag 
coefficient of a cyUndar of radius r spinning about its axis: { « 2ij(wr*0, 
when i| is the viaooaky of the fluid (-0.01-0.1 poise) and t, ram the 
length and radaas, nayectfuly. ForbcUcasc, r»63nm t €»Snm.eothat 
C - 6 • (0.01 [poise]) • (w * 6.5* • 8){nm*} - 6 X HP* pN-nm*. Urns the 
rotary diffusion coefficient ■ k^Tfj - 6 X 10 3 s" 1 . The volume of the 
portal protein Is about twice that of the helicase (r ~ 8 nm, A ~ 10.5 nm), 
so its rotary dUreskm coefficient is about half that of heUcaae. This 
assumes that the viscosity retarding the portal protein's rotation is equiv- 
alent to water, which isprobabry a low estimate.) 



potentials kx^*f = 0, A = r + fr,o> = ir/4) (see Rg. 7): 
w «^ e -^J K + ^j t (2a) 

where x = rB is the coordinate along the DNA charges (the 
diagonal in Pig. 2 c), k is the inverse screening length, and 
b is die distance off-axis (details of die electrostatic calcu- 
lations arc given in the Appendix). Note that in this approx- 
imation the force acting in the translation^ direction is zero. 
The shape of is shown in Fig. 2 a, where the period is 
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The conformations realized by \%B t r) switch among the V, 
values either sequentially or randomly, as will be discussed 
below. 

Fig. 3 b was computed using the method described in 
Appendix B. It shows a typical load torque-angular velocity 
curve for die rotary motion, comparing the perfbtmance 
when the hydrolysis rates at the three sites are independent, 
and when the hydrolysis is correlated so that the flashing 
proceeds sequentially. There is good evidence that, at least 
in tinesin, the hydrolysis rates at different sites are indeed 
correlated (Gilbert et aL, 1995; Hackney, 1994; FesUn and 
Ostar, 1995), so it is quite possible that the sites alternate in 
die rotary motors b e ca us e of steric/elastic coupling. Fig. 3 
shows a comparison between the naecftankal performance 
for three "firing sequences": l-»2-fr3-»l, 3-+2-*l-+3, and 
random. Because of the asymmetry of the flashing potential, 
the 3-+2-+1 sequence outperforms the l-*2-*3 sequence, 
with random firing intermediate between the two. (Note that 
in these simulations precisely a single potential is on at a 
time.) 

Because we have assumed that the angular and axial 
velocities are not independent, the axial velocity along die 
DNA is obtained from the angular velocity times die radius 
by simply multiplying by the sine of the tilt angle of the 
dipoles; because we have chosen the dipole angle as 45°, the 
linear and angular motions are equal. The simulations dem- 
onstrate that the observed linear velocities required for 
helicase progression (-~1<T 2 jim/s) are easily achievable 
with reasonable hydrolysis rates. The portal protein, how- 
ever, requires linear velocities of -1 fim/s, which would 
require hydrolysis rates beyond those observed in other 
ATPases. However, because the portal protein is a much 
larger tredecameric structure, it may have as many as 12 (or 
even more) hydrolysis sites. Wim this many flashing di- 
poles, the model can generate axial speeds in this range. The 
difficulty in simulating the portal protein is the large num- 
ber of possible "firing sequences" with so many hydrolysis 
sites. Given the disparity in performance in the three-site 
helicase model, a wide range of velocities are possible for 
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Hie equations describing thii interaction can be formu- 
lated In terms of the ion concentration at height &G(l0 and 
the coordinate of the movable charge region, y(t). Because 
there are about five or tlx helical tares in traversing a 
biiayer, wc approximate the diffusion of the protons through 
the channel by a set of discrete compartments, in ooe of 
which the interaction with the fixed charges takes place. 
Fig. 4 b shows that the proton flux causes the movable 
charge to oscillate. The reason for this is most easily seen by 
treating the protons as a single compartment adjacent to q. 
This yields a pair of equations similar to the van der Pol 
oscillator. 



Prictional force 

on the 
moveable charge 



Coulombk 
force 



mastic 
force 




(4) 



Proton Association Dissociation 

concentration 
at the fixed 
charge site 

Here t is the time constant for proton relaxation, k is the 
elastic constant of the movable charge, g, r is a "hard-core** 
distance limiting the motion of fi, and (* I0 , ft,, Jtj) arc 
rate constants governing the binding and dissociation of the 
protons to the fixed-charge region. In this approximation we 
have neglected the interaction between the protons. Increas- 
ing the dimensionality of the proton diffusion channel to 
many compartments, or to a continuum diffusion channel, 
does not change the basic dynamics: the proton flux induces 
reciprocating oscillations of the movable charge Q that 
pump the local electrostatic field; in turn, this drives the 
motion of the rotor as analyzed above. 



Reversing the motor's direction 

Whereas helicases are generally unidirectional (either 
5'->3' or 3'-*5'). the portal protein may reverse its direc- 
tion when injecting the phage DN A into the host (Another 
possibility is that die free energy accumulated in p ackaging 
the DNA into the capsid is sufficient to drive injection into 
the host simply by entropic expansion.) The bacterial flagel- 
lar rotor spontaneously reverses its direction of rotation 
stochastically every few seconds. This produces die "runs" 
and "tumbles*' characteristic of bacterial swimming tBerg 
and Hipped 198S); A surprising feature of the model is that 
it can reverse its direction of rotation without changing the 
symmetry of the charge pairs. Simply altering the shape of 



die potentials, for example, by changing the spacing be- 
tween the charged regions, can reverse the direction of 
rotation. - A 
This effect is quite counter-intuitive; it arises because of 

a) the asymmetry of the charge pairs (i*., the else of the 
negative and positive charge regions is not identical) and 

b) the mismatch between the number of stator charge pain 
and the number of charges with which they internet (the 
•Vernier" effect). Prom Bq. 2b we can compute the average 
field experienced by the rotor from all of the charge pairs by 
integrating around the stator circumference. The asymmetry 
of the individual fields and the vernier effect conspire to 
create a net field seen by the assembly whose symmetry can 
switch without switching the individual charge pair symme- 
tries. This is illustrated in Fig. 5 for the case of a stator with 
three charge pain acting on a rotor where two charges are 
within range of their fields. For illustrative purposes we 
have shown die potential as a piecewise linear M wedge M 
whose left and right sides are of length 61 and dr, respec- 
tively (the electrostatic field works as well but Is more 
difficult to visualize). 

Fig. 5 a is a bifurcation diagram showing how altering the 
shape of the potential, but not its symmetry, will reverse the 
direction of rotation. Therefore, to reverse the motor, the 
charge regions do not have to flip their polarity, but need 
only change their separations or angles with respect to the 
rotor. Fig, 5, b and c, shows this effect by illustrating the 
difference in symmetry experienced by a single charge and 
by a pair of charges rotating in the composite field of three 
potentials. 

Because the spacing of the charges is influenced by 
Brownian motion, one can imagine that the operating point 
in Fig. 5 a could stochastically drift between clockwise and 
counterclockwise regions, leading to die stochastic switch- 
ing behavior of the bacterial flagellar motor. We will inves- 
tigate this model for the flagellar rotor more completely in 
another publication. 



DISCUSSION 

We have proposed a model for the mechanochemical oper- 
ation of die rotary protein motors, DNA helicase and 
die portal protein. We have taken our inspiration from the 
recent elucidations of the structures of several rotary protein 
motors. In our view, the common structural features are 
i) the existence of multiple nucleotide hydrolysis sites lo- 
cated circumf erentially and equally spaced in apposition to 
a hydrophobic channel; ii) pairs of charged regions abutting 
die central channel. In at least one of these structures, the 
Fl ATPase, the requisite structures have been identified 
(Abrahams et al., 1994). The model assumes that binding of 
nucleotide to a circumferential site induces a mechanical 
strain in the protein lattice that is transmitted to the charge 
regions. Changing the configuration of the charges with 
respect to one another changes the local electrostatic field 
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let up by the charge regions. Therefore, the negative phos- 
phate charge! that constitute the DNA backbone aee a 
pattern of fluctuating electrostatic fields. These "flashing 
fields" produce a sustained electrostatic torque on the DNA. 
This torque does not require that the fields flash i 
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dally; it arises because each flashing field U asymmetrical. 
To our knowledge, this mechanism for torque generation is 
novel. Pig. 6 summarises the effects producing the rota- 
tional torque. 

Surprisingly, the motor can reverse its direction in re- 
sponse to relatively small alterations in the shape of the 
potential field doe to the charge pair regions. This does not 
require that the fields switch their polarity. Reversal results 
from the "vernier" effect that arises because the number of 
flashing fields is not the same as the number of entrained 
charges. This small shape modification could be controlled, 
for example, by phosphorylation/depbosphorylation (Oold- 
beter and Koshland, 1982). The bacterial flagellar rotor 
reverses stochastically every few seoonds, and this reversal 
mechanism is an attractive candidate for this phenomenon. 

Several authors have invoked the vernier effect in dis- 
cussing the portal protein. In 1978 Hendrix pointed out 
that the pentagonal symmetry of the ioosohedral capsid 
base interfaces the hexagonal symmetry of the co nn e cto r 
(Hendrix, 1978). This led him to suggest that this symmetry 
mi s m a tch might be an essential dement in a rotary motor 
that screwed the DNA into the capsid. He reasoned that the 
mismatch between potential minima on the rotor and stator 
would produce a large number of shallow potentials that 
would offer little resistance to rotation. Although Hendrix 
did not address the mechanism of torque generation, he 
implied that the symmetry mismatch was somehow in- 
volved. However, later meuuranents showed that each 
base pair inserted consumed one-half ATP <Guo et al., 
1987), and the 6-to-5 symmetry mismatch would consume 
too much ATP per revolution (Dube et al, 1993). Recent 
measurements also reveal that the portal protein has a 12- or 
13-fold symmetry, which led Dube et aL (1993) to suggest 
that the rotary motor relies on the di sc rep ancy between the 
tredecameric symmetry of the portal protein and the 10-fold 
helical symmetry of the DNA strand (Dube et al., 1993). 
However, they postulated that insertion was not rotary, but 
resulted from an axial movement In the mechanism of force 
generation we have proposed, the relative motion of the 
rotor and stator has both axial and rotary components. The 
only other model we are aware of that deals with the force 
for p ackagin g is die osmotic pressure model of Scrwer 
(1988). 



FIGURE 5 Reversing the motor, (a) The parameter plane describing the 
shape of the potential. In <+) regions the motor runs clockwise; in (-) 
regions the motor runs co unt e rclo ckwise. The shape of the potential is 
characterized by Al and Ar. (o and c) The shapes of the potential seen by 
each charge on tbe rotor {top) rod the shapes of the potentials teen by the 
assembly of two charges rotating in the field of three potentials {bottom). 
The bottom panels were compnted by inserting the p«««— f«i in the top 
panel into Eq. 2b. In b the symmetry of the single and coe ap o s i te potentials 
are the same, and the rotor moves clockwise. In c the symmetry of the 
single potentials is the same as in n, but the relative steepness of the sides 
has changed. The rotor moves counterclockwise because the co m po si t e 
potential seen by bom charges has symmetry opposite that seen by a single 
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At first glance, the mechanism we have described may 
seem to resemble the "flashing ratchet" mechanism pro- 
posed independently by several authors (Ajdari and Prost, 
1992; Astumian and Bier, 1993; Festdn et al, 1994). How* 
ever, mat mechanism depended on biasing the diffusion of 
the load by flashing an asymmetric potential, and so de- 
pended on Brownian motion in an essential way (i.e., the 
force vanishes as k*T -+ 0). The mechanism we have 
described here generates torque electrostatically, and the 
role of Brownian motion is to act as "lubrication" to jiggle 
the rotor out of local minima (Doering et al., 1994; Mag- 
nasco, 1993). We have also shown that the flashing field can 
also be driven by an ion (e.g., proton) gradient, which 
makes it a candidate for the bacterial flagellar motor. A 
similar mechanism may operate to drive the rotation of the 
y subunit of ATP synthase (Abrahams et al., 1994). 

How can the model be tested? The model predicts a 
specific form for the load-velocity relationship that appears 
to be different from those measured for other molecular 
motors (Berg and Turner, 1994; Svoboda and Block, 1994). 
Moreover, the statistical properties of the motor* s morion 
are distinct from other mechanisms (Samuel and Berg, 
1995; Svoboda et al., 1994). The fact that the motor is 
electrostatic suggests a more direct measurement: an alter- 
nating electric field imposed across the motor would inter- 
fere with torque generation when its frequency is near that 
of the flashing field. Although quantitative load-velocity 
and variance curves remain to be measured, the model 
makes definite predictions as to the quantities mat need to 
measured, and provides a novel mechanism for reversible 
rotary motors. 

APPENDIX 

A; The motion of a helical array of charges 
driven by a flashing dipole 

We consider • cylinder with reditu R on which an amy of dipoies 
centered et c 0 ere positioned with angle e) with respect to the tangent 
to the cylinder (see Fig. 7). Thus, e) » 0 are dipoies along the radios of 
the cylinder with no z co mp onent. Let d be the half-distance between 



the two charges; then simple tri tonometry shows that the position of tha 
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where r Is the diameter of the DNA hefix (r < K). o is the pitch, and the 
variable 6 denotes the rotational oo mp one oi of the helix. Let the position 
of the charges along the helix be given by f • c + /Cq» where c hi between 
0 and to sndj is an integer. Thus, the precise position of die helix at any 
given time is specified by (0. t). Because t - 0 is Identified with t - t» 
die coordinates of the helix are descri b ed by a pair of values on the torus; 
hence the potential felt by the charges on die DNA is a doubly periodic 
ftiocftoOi To compntp thii ftutctioo wo boGb (hst the potBClfaJ nfels due to a 
s dipole at posi t ion jj is given by 
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(Jj.yj.ii) (A4) 
= (rcos afc + jzj - d, r sin a(e +/Zo) - 0, t +/zo). 
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FIGURE 7 Geocnetryof the dipok locations with napect to the negative 
charges oo the MIA backbone. 



wherein, - 1 If the J1tk dipofo is M on f ** otherwise it is zero. Here £f(0 are 
independent white noise variables. To simplify things, we assume as, Is 
zero except for exacfyone value of I and that the rate of tnnthloo from 
one state to the next fi^onlermmod by r. u te 

we assume no scrojosjj and take a - X b • 8, c -4U, r - 20, c - 0 J, 
- 6, - 5. The is is agraUon was carried out using a fixed stop size of 
0.01 for 200j000 steps using the program xpp. 



B. Approximations to the flashing field equations 

We shall modal me system as a periodic sequence of charges (DNA 
b a ckbone ptjospbeses) moving through a p e ri odic array of **off-«jris- 
saeenc<Mipok , * potentials, each governed by Eq. 2. The phase of the 
underlying potential Is allowed to change because of ATP-geaented coo- 
forniational changes moving or "turning on and ofT the dipoka. The 
rwtentiaJduetoeschdspoteuiti^^ 
DNA is 

V^=Sv^-2nrj). (Bl) 
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where (is the friction coefficient, V%x> I) It the time oopandent potential, 
"MAtfoy to*** to «W>er V,C«X or V,C»), with swls^ rukss to 
be speci fi ed below, sod x(r) Is a Osussta a white noise process satisfying 

<«»)>- 0, (Mm -*>-*). (B4) 

Pot n o t a ti o n s! si mp l i cicy wa divide the stochastic differential equation 
mro ng h by C* identify the diffusion coe ffi cien t D m fcgT/{, and resells dm 
p ote nt ia ls by so the equation becomes 



- 3? --V(X,l)+^5«l). (BS) 

The switching dynamics am inodesod as random, Marfcovian ttanshsons 
between die different oo nfl gu ra tions. This means that the watting time In 
each ww fl guration is a P oi sso rt distributed rsmtflm variable, independent of 
all the other dynamics. The time acsJeof meswisdiM todiemvemof dm 
switching rase, which win be denoted a; Le*. a Is the rate of jumping 
between sub s equen t configurations and a' 1 Is the average time between 
Bwncbos. 

These stoch a sti c dynamics are then fully described by a set of coupled 
Fokker-Plsnck equatkms tor die evolutira 

of finding the poomtial in configuration n and the charge army at position 
equations have (he form 
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/F,-ot Via Via \/P,<JM)\ 
= Via Ft- a Via HlfeOl 
\ Via Via Fj-a/W,')/ 

where a is the rate of switching and die differential op erat or F m is 

<m this notation, an open derivative operator like Mdx operates on every- 
thing to its right) For the case of sequential flashing in order l-»2-*3-~»l 
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a constant value as the 
velocity is simply 
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This formulation It for the "no-toad" situation. If m external load torque, 
/, is applied to the DNA, tea the equivalent particle dyntmict b the tune 
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where K + 1 it to be interpreted u 1, and 0 U to be Interpreted at K. 
In the absence of switching, this discretization supports a troe equilib- 
rium stationary state with vanishing currents tod Boltzmann probability 
distributions ■ 

The numerical eolations to these equations were canied out using 
MatLab on a Sun Span-station 20. The code is available via internet ftp 
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